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Abstract. We present a new model for soft gamma-ray repeaters based on a quark star born with temperatures 
above the critical value (T c ) for the onset of the colour-flavor locked superconductivity. The quark star then 
quickly cools below T c , expelling a fraction of the surface magnetic field via the Meissner effect. We show that if 
a small fraction (< 10%) of the surface magnetic field (10 14 — 10 15 G) is expelled, it quickly decays via magnetic 
reconnection and heats up the quark star surface to temperatures > 10 9 K. Created (e + ,e~) pairs annihilate 
into gamma rays emitted in a giant burst (the first burst in our model), with a luminosity of ~ 10 45 ergss _1 . 
Subsequent bursts result from the restructuring of the surface magnetic field following the formation and relaxation 
of a vortex lattice which confines the internal magnetic field. During this phase, energy is sporadically released 
as a consequence of magnetic reconnection events in the entangled surface magnetic field as it evolves into a 
smooth, more stable, configuration. The star eventually enters a quiescent phase in which energy is continuously 
supplied by vortex annihilation at the surface. As the star spins down, the outermost vortex lines will be pushed 
to the surface where they annihilate and release their confined magnetic field. We show that the corresponding 
luminosity is L v ~ 10 36 ergss -1 for a typical soft gamma-ray repeater spinning with a period of 8 s and a surface 
magnetic field not exceeding 10 15 G. Our model can be applied to any situation where a T > T c quark star is 
generated. We discuss the connection between anomalous X-ray pulsars and soft gamma-ray repeaters in the 
context of our model. 
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1. Introduction 

Soft 7-ray repeaters (SGRs) are sources of recurrent, short 
(t ~ 0.1s), intense (L ~ 10 3 — 10 4 I/Edd) bursts of 7-ray 
emission with a soft energy spectrum. The normal pattern 
of SGR activity are intense activity periods which can last 
weeks or months, separated by quiescent phases lasting 
years or decades. The five known SGRs are located in 
our Galaxy, or, in the case of SGR 0526-66, in the Large 
Magellanic Cloud. The two most intense SGR bursts ever 
recorded were the 5 March 1979 giant flare of SGR 0526- 
66 (Mazets et al. 1979) and the similar 27 August 1998 
giant flare of SGR 1900+14. The peak luminosities of these 
events ( ~ 10 6 — 10 7 LEdd) exceeded the peak luminosities 
of "normal" SGR bursts by a factor > 10 3 . In Table^we 
summarize the burst properties of these giant flares. 

Several SGRs have been found to be X-ray pulsars with 
an unusually high spin down rate of PjP ~ 10~ 10 s _1 , 
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usually attributed to magnetic braking caused by a super- 
strong magnetic field B > 10 14 G, which implies that 
SGRs are magnetars (Golenetskij et al. 1979; Duncan & 
Thompson 1992, Kouveliotou et al. 1998, Kouveliotou et 
al. 1999). In the magnetar model, the magnetic field is 
the likely provider of the burst energy, since it is the 
dominant source of free energy in the star. A common 
scenario is that stresses build up in the magnetic field 
and create a quake in the crust of the neutron star which 
ejects hot plasma Alfven waves through its rigid magneto- 
sphere (Thompson & Duncan 1995; 1996). The magnetic 
field of such a star would have grown to magnetar-scale 
strengths because of strong convection during the collapse 
of the proto-neutron star core (Duncan & Thompson 1992; 
Thompson & Duncan 1993). 

In this paper, we propose an alternative model where 
the SGR activity is produced by phenomena occurring in 
a hot (T > T c ) quark star. In §|2we discuss the cooling 
timescales of the quark star and the onset of colour su- 
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Table 1. Giant burst properties. 



Object 


SGR 0526-66 


SGR 1900+14 


Active periods 


1979-8J 


1979, 1992, 1998-99 


Giant Burst 


5 March 1979 


27 August 1998 


Precursor: 






Duration, s 


no data 


~ 0.05 


kT (keV) 




~ 20 


Hard 7-ray spike: 






Duration, s 


~ 0.25 


~ 0.35 


Peak luminosity, ergs s~ 


1.6 x 10 45 


^ 3.7 x 10 44 


Energy release, ergs 


1.3 x 10 44 


^ 6.8 x 10 43 


kT (keV) 


~ 246 


~ 240 


Bright X-ray emission: 






Duration, s 


~ 180 


~ 370 


Pulsation period, s 


8.1 


5.16 


Energy release, ergs 


3.6 x 10 44 


5.2 x 10 43 


kT (keV) 


~ 30 


evolves, 31.2 -> 28.9 



Note 1. Data sources: Cline et al. 1980; Fenimore, Klebesadel, 
& Laros 1996; Feroci et al. 2001; Hurley et al. 1999a; Ibrahim 
et al. 2001; Mazets et al. 1999. "No data" means that Mazets 
et al. 1979 have not plotted any data for the last second before 
the burst (the time resolution for their observations of the pre- 
burst background is apparently > Is). 



where a c = g 2 /An is the QCD fine structure constant, 
g is the quark-gluon coupling constant, Y e = n e /rib is 
the ratio of electron and baryon numbers, no ~ 1.7 x 
10 38 cm -3 is the nuclear matter saturation density, Tg is 
the temperature in units of 10 9 K, and taking the specific 
heat per volume as (Iwamoto 1982) 



C q = 2.5 x 10 20 (nb/no) 2/3 r 9 ergscm^KT 1 , 



(2) 



we can estimate the timescale for neutrino cooling from an 
initial temperature Tg^ to a final temperature Tgj <C Tgy. 



Tcool ~3X lO 2 ^ 1 ^ 1 / 3 I — 



-1/3 



Ta 



(3) 



According to this estimate a critical temperature of, say, 
T c ~ 10 11 K will be reached within r coo i ~ 0.4 ms. Thus, 
the star will rapidly undergo a phase transition to a colour 
superconducting state (except in the cores of the vortices 
formed in response to the rotation of the star; see § 0} . 



perconductivity. In § we show how giant flares such as 
the 5 March 1979 event could be powered by the expul- 
sion of a fraction of the magnetic field due to the Meissner 
effect after the onset of superconductivity. We also calcu- 
late the light curve and temperature evolution for such a 
burst in our model. The remaining magnetic field in the 
star is confined to vortices, which will evolve into a lat- 
tice configuration, as discussed in §0] Subsequent weaker 
SGR bursts may be due to the reorganization of the ex- 
ternal magnetic field, and the luminosity of SGRs in their 
quiescent phase may be mostly powered by the annihila- 
tion of vortices at the stellar surface, as detailed in § [5] 
After decades of quiescence, a crust may form, leading 
to occasional crustquakes which could power later burst 
events, and additional events could be due to rare im- 
pacts; such events are discussed in § In § we briefly 
discuss known SGRs within our model, and in §[S]we dis- 
cuss the suggested connection SGRs have with anomalous 
X-ray pulsars (AXPs) in relation to our model. Finally, 
we conclude in § 

2. Cooling timescale and onset of colour 
superconductivity: the CFL star 

Present estimates of the critical temperature for onset 
of colour-flavor locked (CFL) superconductivity suggest 
T c ~ 10 12 K (Rajagopal & Wilczek 2000) whereas early 
estimates indicated T c ~ 10 9 K (Bailin & Love 1984). The 
quark star we take to be born with an effective temper- 
ature higher than T c . Once formed, the quark star will 
cool rapidly through neutrino emission. Using the result 
of Haensel (1991) for the neutrino emissivity, 



£„ = 2.2 x 



io 2 vn 1/3 K/rco)T 9 e 



ergs cm s 



(1) 



3. Giant Burst: Expulsion and decay of the 
magnetic field 

In an ordinary superconductor there is a thermodynami- 
cal critical field H c (or H c i for type II superconductors), 
determined by the free-energy difference between the nor- 
mal and superconducting states in zero field, below which 
a magnetic field will be screened from the interior of the 
superconductor: this is the Meissner effect. The existence 
of the magnetic Meissner effect in colour superconductors 
is a matter of some debate (Alford, Berges & Rajagopal 
2000, Sedrakian et al. 2001, Iida & Baym 2002) but we will 
proceed on the assumption that there is a Meissner effect 
in the CFL phase, and that at least part of the magnetic 
field is expelled from the surface layers of the star. 

The model we consider assumes that the star 
intially consists of uniform density strange quark 
matter with an associated surface electric field 
|Alcock, Farhi, fc Olinto 1986| In this situation the 
calculations of |Page &: Usov (2002) show that the surface 
of the star cools faster than the core because of thermal 
emission of photons and e + e~-pairs from the surface. 
Therefore the CFL phase transition occurs first in a thin 
layer of thickness S -C R (R is the radius of the star) close 
to the surface. The magnetic field decay discussed later 
will thus heat the surface leading to thermal emission of 
photons. 

The penetration depth A for the magnetic field is of 
the order 1 fm (see Iida and Baym 2002) so the magnetic 
field is negligible throughout the superconducting layer 
provided that S ^> A ~ fm. Note that the thermodynami- 
cal critical field is H c ~ 10 19 G (Iida & Baym 2002), much 
larger than the field strengths of interest to us. 
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3.1. Expulsion timescale 

The time scale for the expulsion is determined by the 
conductivity of the normal, non-superconducting state ac- 
cording to 

Ana e \d 2 ( B \ 2 



' cxp 



where a c \ is the electrical conductivity in the normal state 
while B, c and d are the magnetic field strength, the 
speed of light and the thickness of the superconducting 
layer, respectively. For quarks, this quantity was found by 
Heiselberg & Pethick (1993) to be given by 

^-^^'(aoow)' 7 ^ 101 " 8 " 1 ' (5) 

where fi is the quark chemical potential. For a s = 0.1 
(strong coupling constant), a c \ varies from ~ 10 17 s _1 



Quark Star with T > Tc [u,d,s] 




Quark Star with T < Tc (CFL) 




Vortices 

Fig. 1. Schematic view of the initial and final stages in the 
development of the SGR in our model. The QS, born with 
T > T c quickly cools into a CFL star. In the early stages 
of this transition, a fraction of the surface magnetic field is 
expelled by the Meissner effect, inducing magnetic recon- 
nection events leading to the main burst. In the late stage, 
once the vortices form (parallel to the rotation axis), the 
system evolves into a configuration where most of the mag- 
netic field is confined to the vortex lines. Magnetic recon- 
nection events between the complex surface magnetic field 
and the external part of the magnetic field coupled to the 
vortex lines lead to occasional flares (subsequent bursts). 
In the SGR quiescent phase, magnetic energy is released 
as vortex lines are being continuously pushed to the sur- 
face and annihilate there. In the more general case where 
the magnetic field is not aligned with the rotation axis, 
we expect an even more complex behaviour/evolution of 
the field leading to more reconnection events. 



for T = 10 11 K to 
electron contribution to the electrical conductivity can be 
estimated from the classical result 



^el, electron 



m e v e 



(6) 



where a ~ 1/137 is the QED fine structure constant, n e 
is the electron number density (m e is the electron mass), 
£ e is the electron mean free path, v e their mean thermal 
velocity, and taking the ultrarelativistic limit where v e ~ c 



and m e is replaced by T, using £ e 
estimate, this gives 



-1/3 



as a crude 



cr cl ~ 7.1 x lO- 



ii ( Y e 



V 10- 



2/3 / \ 2/3 / w A/r x r x 

f n h \ (I MeV\ _ x 

s ■ ( 7 ) 

n 0/ 



— 



More refined estimates (Baym & Heiselberg 1997; Arnold, 

Moore & Yaffe 2000; Shovkovy & Ellis 2003) give cr i ~ 
10 23_ 1Q 24 s -i t With R = 1Q i9 q we btain 



' cxp 



10" 



/ Q"cl 

VlO 17 s" 1 



B 



10 3 cmy \10 15 G 
is r oxp of th 

onds or less since cV ax ~ R ~ 10 km 



(8) 



so a reasonable estimate is r exp of the order of a few sec- 



3.2. Energetics 

The expelled magnetic field would quickly lead to mag- 
netic reconnection in the magnetosphere. Since reconnec- 
tion typically occurs at a fraction of the Alfven velocity, 
the growth time of the instability can be estimated to be 
~ 10~ 5 — 10~ 4 s which is comparable to the ~ 2 x 10~ 4 
s rise time of the March 5 event (Mazets et al. 1979; 
Paczyhski 1992). The long expulsion time (see Eq. [S] for 
$max. ~ R) implies that the reconnection events contin- 
ued for many times the Alfven time scale which we suggest 
could account for the longer duration of the hard transient 
phase of the burst. 

The energy released is 



E 



M 



1.7 x 10 47 ergs 



(3B 
10 15 G 



Rqs 
10 km 



(9) 



where f3 is the fraction of the surface magnetic field ex- 
pelled with the subsequent reconnection events (the sub- 
script "M" stands for main) . Locally we expect the energy 
from the magnetic reconnection event to be rapidly con- 
verted to a thermal photon-pair plasma (e.g. Thompson 
and Duncan 1995, section 3) which heats up the surface. 
The cooling is defined by the rate at which the photons 
escape this heated region. However, since these reconnec- 
tion events occur very close to the surface of the star, we 
expect the photons to be trapped, and the cooling time 
to be long enough to allow for thermalization. The corre- 
sponding thermal temperature is 



T M - 4.0 x 10 9 K 



(3B 
10 15 G 



1/2 



10 km 



1/4 



(10) 
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Fig. 2. The light curve expected in our model. The four 
curves (from top to bottom) correspond to /3 = 0.5, 0.2, 
0.1, and 0.05. 



which is the maximum surface temperature that can be 
reached. Even smaller surface temperatures (when cooling 
is faster than heating) are enough to trigger the mecha- 
nism of thermal emission, as we discuss next. 

3.3. Light curve 

Usov (2001) showed that creation of e + e~ pairs by the 
Coulomb barrier at the quark star surface is the main 
mechanism of thermal emission from their surface at the 
temperature Ts < 5 x 10 10 K. Created e + e~ pairs mostly 
annihilate in the vicinity of the quark star into 7-rays, and 
Usov (2001) argued that the light curves of the March 5 
1979 and August 27 1998 events may be explained in a 
model where the burst radiation is produced by the bare 
surfaces of such stars heated up to ^ 2 x 10 9 K by impacts 
of massive cometlike objects. He also points out that any 
other mechanism which quickly heats up the surface can 
explain the events, so Eq. I|10|) indicates that the Meissner 
effect and subsequent decay of the expelled magnetic field 
can power a burst. To give further evidence for this we 
have computed the light curves in our model. 

We consider a thin layer of thickness S <C R close to 
the surface of the star which is heated by the release of 
magnetic energy E^ec from the decaying magnetic field. 
The temperature in the layer (taken to be isothermal) is 
governed by the equation 



VgC— = Q- V 5 e v , 
at 



(11) 



where V$ = 4irR 2 5 is the volume of the layer, and Q for 
times < t < Td oc is given by the energy released per unit 
time by the decay of the magnetic field, Q ~ E& ec /T& cc , 
while for times t > Td cc it is given by the luminosity of 



m e c 2 + kT is the mean energy of the electron-positron 
pairs created, 



/± = 1.6 x 10 39 T 9 3 exp 



11.9 



J(C) cm^s" 1 



is the flux of pairs per unit surface area, 



J(0 = 



c 4 



lC 3 ln(l + 2C" 1 ) 7T 5 
3 (1 + 0.074C) 3 + IT (13.9 + C) 4 ' 



(12) 



(13) 



and C = 20/Tg. The specific heat per unit volume C is 
the sum of the contribution from the electrons (Blaschke, 
Grigorian, & Voskresensky 2001) 



C e ~ 5.7 x W 19 Y 2/3 {n b /n Q ) 2/3 T g ergs cm" 3 K~ 



(14) 



and the heat capacity of the quarks in the CFL phase is 
modified from Eq. J5J) according to 



C q = C q f(T/T c ), 
where C q is given by Eq. J2J and 
1 



ft \ 3 - 2 
t [x = — exp , 

x M 0.57a; 



(2.5- 1.5a; + 3.6a; 2 ), 



(15) 



(16) 



see Horvath, Benvenuto, & Vucetich (1991). Also, the neu- 
trino luminosity is suppressed in the CFL phase by the 
Boltzmann factor exp(— A/k&T), where A is the energy 
gap, i.e. 



e v = € u g(T/T c ), 
where 



g{x) = exp 



0.57a; 



(17) 



(18) 



the e + e pair emission, Q 



-4ttR ! e±/±, where e± 



and we have assumed that the temperature dependence 
of the energy gap is A(T) = A(T = 0)(1 - T/Tc) 1 ' 2 (e.g. 
Carter & Reddy 2000), and used the standard relation 
k B T c ~ 0.57A(T = 0). 

We have solved equation (|ll|l numerically with the ini- 
tial condition T = 10 s K. For the parameters we used the 
values a c — 0.1, rib — 5n , Y e = 10~ 4 , R = 10 6 cm, 
5 = 10 3 cm, B = 10 15 G, r doc = 0.1. We varied the frac- 
tion of the magnetic field expelled by the partial Meissner 
effect, using (3 = 0.5, 0.2, 0.1 and 0.05 as representative 
values. The light curves are shown in Figure [2] The maxi- 
mum luminosity for (3 — 0.1 is 5 x 10 45 ergss -1 , so the de- 
cay of 10 % of a surface magnetic field of 10 15 G is enough 
to power the main burst. The light curve and time evolu- 
tion of the surface temperature for this case are shown in 
Figure [3j 

We note that since the CFL phase transition occurs 
only once, there will be only one giant burst in our model. 
We also note that our model is too simplistic to repro- 
duce the periodic pattern (e.g. the 8.0 s period in the 
March 5 event) overimposed to a smooth exponential de- 
cay. However, given the partial Meissner effect and the 
resulting random reconnection fronts we expect isolated 
spots on the surface of the star to be heated and release 
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4. Subsequent bursts: Vortex lattice dynamics 

4.1. Vortex formation and relaxation 

The CFL front quickly expands to the entire star followed 
by the formation of rotationally induced vortices, anal- 
ogously to rotating superfluid He 3 (the vortex lines are 
parallel to the rotation axis; Tilley&Tilley 1990). Via the 
Meissner effect, the magnetic field is partially screened 
from the regions outside the vortex cores. The system now 
consists of alternating regions of superconducting material 
with a screened magnetic field and the vortices where most 
of the magnetic field resides. 

The vortex relaxation time, which is within an order 
of magnitude similar to the random diffusion timescale, 
can be estimated to be 



1.5 



1.0 



0.5 



0.0 

-200 200 400 600 800 

t(s) 

Fig. 3. Upper panel: The light curve for j3 = 0.1. Lower 
panel: The corresponding time evolution of the surface 
temperature. 



energy as described above. These hot spots would pulse at 
the same rate as the rotation of the star, and thus, would 
appear periodic. We further note that a few of these hot 
spots could appear simultaneously at random locations on 
the surface of the star. This means that each of the pulses 
(superimposed on the smooth curve) could consist of even 
smaller subpulses. Multiple small hot spots should pro- 
duce many harmonics in the lightcurve. In other words, 
if our model is a correct representation of SGRs, observa- 
tions could constrain the number of hot spots. 

Finally, we note that given the rapid cooling time scale 
of a quark star, as discussed in section 2, the CFL phase 
transition should happen immediately after the formation 
of the quark star. However, this does not necessarily mean 
that the giant burst will follow immediately after the core 
collapse in a supernova since the quark star can in princi- 
ple be formed also much later. 



2?1 

A v c 



(19) 



where A v = (nR 2 /N^) 1 / 2 is the vortex mean free path, 
and iV v is the number of vortices in the entire star, 



<U - 10" [ ^ 



M/3 



300MeV 



(— V 

V 10km / 



(20) 



(Iida & Baym 2002), where /j, is the quark chemical poten- 
tial related to the density. Using the above values, we thus 
estimate a relaxation timescale of ~ 11 min for P = 5 s. 
Note that this is within an order of magnitude consistent 
with the relation between lattice formation time and rota- 
tion period found in numerical simulations of vortices in 
rotating Bose-Einstein condensates (Tsubota, Kasamatsu, 
& Ueda 2002). 



4.2. Magnetic field restructu ration and reorganization 

Given the complicated structure of the resulting surface 
magnetic field, following relaxation, the latter may then 
suffer from frequent magnetic reconnections and thus ac- 
count for the subsequent bursts. This phase could last 
for days or months, depending on the magnetic field 
strength, the period and the fraction of the field which 
remained entangled in the period immediately following 
the giant burst. These random reconnection events would 
bear many (temporal and spectral) similarities to the 
main burst but we expect them to be less energetic as 
the magnetic field slowly decays and weakens. Eventually, 
the magnetic field evolves into a stable configuration (see 
Figure HJ after which the star enters a quiescent phase. 
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5. Quiescent phase: Vortex annihilation 

The number of vortices decreases with the spin-down of 
the star, and can be formulated as (extrapolating from 
results established for cylindrically rotating superfluids): 1 



dN v 



No 



R 



No 

n 



(21) 



As the star spins down, the outermost vortex line will be 
pushed to the surface and annihilate there (Ruutu et al. 
1997). The corresponding luminosity is 



dE v 
dt 



dN v 
' dt 



dN v dn n 

= e— —— = e v N a 



dfl dt 







(22) 



where e v = (q + e m ~ e m , is the energy density per vortex 
which consists of the rest mass energy (eq) and the con- 
fined magnetic energy (e m ). This means e v Nq ~ E mas — 
1.7 x 10 47 ergs (B/10 15 G) 2 (i?/10km) 3 (recall that most of 
the magnetic field/energy resides in the vortex and that 
e m >> e v ). The star continues to loose rotational energy 
due to the electromagnetic radiation losses which allow us 
to write in the simplest approximation, 



n 



4B 2 R 6 n 2 
9ic3 



B 



10 15 G 



3.24 x 10 

R ^ 



^V 1 x 



10k 



m 



Al 



© 



Al 



(23) 



where / is the star's moment of inertia (Manchester & 
Taylor 1977). 

Combining equation H22|) and equation l|23|) . the lumi- 
nosity induced by vortex annihilation is 



L v = 5.51 x 10 37 ergss 1 x 

— Vf— Vf- 

10 15 G/ \l0km.) \P 



Al 







Al 



(24) 



This can be compared to the luminosity from the dipolc 
radiation. 



Ldip — ~ 



4B 2 R 6 n 4 



2.6 x 10 36 ergss _1 x 



9c 3 

10 15 G ) I 10km ) \p 



6 / -, \ 4 

Is 



(25) 



We note that for P > 1 s we can still account for the 
observed SGR quiescent luminosity with B < 10 15 G while 



1 In rotating superfluids the density of vortex lines in equi- 
librium rotation is n v = 2£^, where k = vh/2m is the cir- 
culation of a vortex with quantization number v for atoms 
of mass m (Ruutu et al. 1997); h is the planck's constant. 
In the continuum limit a totally filled cylindrical container 
with radius R would have iVo = irR 2 n v lines. Interactions 
with the lateral walls give rise to an annular vortex-free re- 
gion along the wall. Its width A is of order the intervortex 
distance r v = (K/2nD.) 1 ^ 2 (expressed here as the radius of the 
Wigner-Seitz unit cell of the vortex lattice). As a result, the 
total line number JV(fl) ~ No(l — 2\ V /R) is always less than 

ivb(n). 



10 15 G (our model) 



,SGR0525-66 

*SGR1806-20(2) 
*SGR1 806-20(1) 
„SGR1 627-41 



t SGR1900+14 



10 15 G (dipole radiation) 




10 



P [s] 



Fig. 4. The X-ray luminosity versus period for the four 
SGRs for which these values have been measured (from 
Kaplan 2000, and references therein). The solid lines are 
predictions given by equation (|24l) for stars with magnetic 
fields of 10 15 G (upper curve) and 10 14 G, for R = 10 km 
and M = IMq. The dashed lines are predictions given 
by equation j^jl f° r stars with magnetic fields of 10 15 G 
(upper curve) and 10 14 G, for R = 10 km. 



much higher fields are required to power the luminosity 
with dipole radiation (Figure 4). 

In our model, the quiescent phase emission consists 
of two sources of energy, related to vortex annihilation 
and dipole radiation. It would be interesting to compute 
the corresponding spectra and compare them to the now 
well-established best fit model of SGR spectra which is a 
composite of blackbody and power-law emission (Kaplan 
2002; Kulkarni et al. 2003). This is beyond the scope of 
this paper and is left as an avenue for future work. 

6. Later bursts: Crustquakes and debris impacts 

6.1. Crust formation and crustquakes 

A crust of hadronic matter suspended above the quark 
star surface (Alcock, Farhi, & Olinto 1986) might even- 
tually form, leading to possible starquakes. Since known 
SGRs are solitary objects, such a crust should have formed 
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due to fall-back of supernova ejecta. Also in the Quark- 
Nova (Ouyed et al. 2002) picture, where a hot quark star 
is formed, high accretion rates are expected from the 
fallback material (Keranen&Ouyed 2003). Later bursts 
could be driven by energy released during crust frac- 
tures that result from magnetic stresses as described in 
Thompson&Duncan (2001). In the case of the neutron star 
crust, the corresponding energy is 

x (f^) 1 ' 6 f^^Vergs, (26) 

where ip is the shear strain, p^ the nuclear saturation den- 
sity and AR^ is the extent of the fracture in the crust. If a 
similar phenomenon can be expected for quark star crusts, 
we expect later (crust induced) bursts with energies as 
high as ~ 10 40 ergs (the maximum density for a quark star 
crust is close to the neutron drip value, pemst ~ 4 x 10 11 
g cm -3 see e.g. Glendenning 1997). 

6.2. Hard emission bursts 

Rare (< 1%) subsequent bursts of hard emission (Woods 
et al. 1999) could be due to random debris impacts. 
Again, as an example, such debris could naturally re- 
sult during disk (and the subsequent planetary) formation 
around the quark star following an expulsion of the crust 
(Keranen&Ouyed 2003). 

7. SGRs in our model 

7.1. SGR 0526-66 

Applying our model to the SGR responsible for the famous 
1979 March 5 giant burst, and for P = 8 s, we find the lu- 
minosity in the quiescent phase to be 8.6 x 10 35 ergss _1 
for B — 10 15 G. The corresponding minimum age, making 
use of equation 1)23(1 . is ~ 3kyr. Hence, our model does 
not require extreme magnetic fields (> 10 15 G) and pre- 
dicts ages of at least a few kyr, in accordance with age 
estimates of the nearby SNR N49 (e.g., Vancura et al. 
1992). However, we note that the association with SNR 
N49 is tenuous, and the large offset between the two ob- 
jects would still be a challenge if the CFL phase transition 
occurred as recently as 1979 (see H8.4JI . Following its long 
quiescent phases SGR 0526-66 could have acquired a thin 
crust (SJjjjJ which could explain why it currently shows X- 
ray characteristics similar to AXPs (Kulkarni et al. 2003), 
two decades after becoming quiescent. We speculate that 
the 1979 March 5 giant burst, detailed in Tabled could 
have been a signature of a CFL phase transition and as- 
sociated physical processes, as described in This SGR, 
displaying a giant burst and a subsequent active phase 
followed by a long quiescence, fits best within the picture 
outlined in this paper. 



7.2. SGR 1900+14 

Because of the presence of SGR burst activity prior to 
the main event (Hurley et al. 1999), the 1998 August 27 
event in SGR 1900+14 would be best explained by, e.g., 
restructuring of the magnetic field (§ 0J, crustquakes or 
debris impact (§0). 

8. Discussion 

8.1. The suggested SGR-AXP connection 

It has been debated whether anomalous X-ray pulsars 
(AXPs) and SGRs are connected, with both classes of 
systems being magnetars (e.g., Chatterjee, Hernquist & 
Narayan 2000; Gavriil, Kaspi, & Woods 2002). AXPs dis- 
play persistent strong X-ray emission with P ~ 6 — 12 s 
pulsations. The term "anomalous" comes from the fact 
that the X-ray emission is not powered by rotational en- 
ergy or by accretion from a companion star (AXPs are 
solitary objects). The likely association of three AXPs 
with supernova remnants indicate that they are young 
(t < 10 yr) systems with unusually fast spin-down rates 
(Gaensler et al. 2001). Like SGRs, the AXPs are rare ob- 
jects (five confirmed cases are currently known), and they 
have a similardistribution of rotational periods. 

If AXPs are magnetars, their emission is most likely 
powered by the decaying magnetic field. Recent observa- 
tional results have suggested a link between the two classes 
of objects: Gavriil et al. (2002) report SGR-like X-ray 
bursts from the AXPs IE 1048.1-5937 and IE 2259.1+586, 
and Kaspi et al. (2003) report a major SGR-like X-ray out- 
burst from IE 2259.1+586. Furthermore, Kulkarni et al. 
(2003) report Chandra observations of SGR 0526-66 in a 
quiescent phase showing that the object has X-ray prop- 
erties similar to an AXP. These similarities have been sug- 
gested to favour a common magnetar model for AXPs and 
SGRs (see Kaspi 2004 for a recent discussion). However, it 
remains to be seen if these two classes of objects are simi- 
lar physical systems that could perhaps be linked through 
a simple evolutionary model, or whether they are in fact 
disparate. 

In our model, SGRs are strange matter quark stars 
that have undergone the phase transition into a colour 
superconducting state. Based on this, we cannot identify 
a simple connection between AXPs and SGRs, but our 
model merely suggests that they are not closely connected. 
We discuss below how missing glitch activity and inferred 
transverse velocities of SGRs could hint in favour of our 
picture. 

8.2. Glitches in our model 

Glitch activity is apparent in AXPs, but seems to be ab- 
sent in SGRs, see Woods et al. (2003a and b) and Kaspi 
et al. (2003). Vortices are a natural outcome of the CFL 
phase transition, and the interaction between the vortices 
and the crust could in principle lead to glitch activity (see 
Alpar 1991). If only a tiny crust is expected for a CFL star, 
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glitches would be unlikely unless a glitch mechanism which 
does not involve the crust-core interaction is involved. A 
possibility for such a mechanism has been suggested by 
Alford, Bowers and Rajagopal (2001). They propose that 
in a certain range of densities quarks may form Cooper 
pairs with nonzero momentum, leading to pairing energy 
gaps which vary periodically in a crystalline pattern, and 
this structure opens up the possibility of vortex pinning. 
However, the details of this mechanism have not yet been 
worked out. Thus we do not expect glitch activity to oc- 
cur in SGRs in our model except maybe in a case where 
a sizeable crust has been formed. This might be possible 
e.g. if the hadron-quark phase transition has happened 
very quickly after the supernova explosion and the crust 
has been formed from the possible fall-back matter of the 
supernova ejecta. 

8.3. Transverse velocities 

At least half of the AXPs are now confirmed to be lo- 
cated near the centers of supernova remnants (SNRs) (see 
discussion by Gaenslcr 2002). This infers transverse veloc- 
ities of < 500kms _1 and an age of < lOkyr. Furthermore, 
three out of the five confirmed SGRs are located near 
SNRs, but the SGR separation from the SNR centers 
and the SNR ages imply transverse velocities of order 
> 1000 km s -1 , which is higher than what would be ex- 
pected from a SN explosion alone. In our model, an neu- 
tron star undergoes a phase transition into quark matter. 
This phase transition is accompanied by a neutrino emis- 
sion (see e.g. Iwamoto 1980). If the neutrino emission is 
not fully isotropic, this may lead into a high transverse 
velocity of the quark star, explaining the large offsets be- 
tween some SGRs and apparently associated nearby young 
SNRs. 

8.4. Formation scenarios 

The transition from hadronic matter into quark matter 
can in principle happen immediately during or after the 
supernova explosion, but also very much later than that. If 
formed quickly after the supernova occurred, the emitted 
gamma rays would be absorbed by the expanding super- 
nova ejecta, and there would be no detectable gamma- 
ray signature of the transition from neutron star matter 
into quark matter. In this case, the resulting quark star 
would be born in a relatively dense evironment and might 
form a crust from accreted material. We would then expect 
subsequent SGR-type bursting behaviour to be caused by 
mechanisms discussed in § 6. As discussed above, because 
of the crust, the star could show glitch activity and could 
in this respect have characteristics more similar to AXPs 
than other SGRs. 

The hadron-quark transition may also happen much 
later than the supernova explosion. Quark star formation 
can, for example, happen as a result of accretion or spin- 
down of a neutron star. In this case, the CFL phase transi- 



tion should be accompanied by an observable gamma ray 
emission (as a giant burst in our picture), as we specu- 
late is the case for SGR0526-66. In this second scenario, 
the quark star would remain bare or have only a thin 
crust due to low accretion (since the objects are solitary). 
Therefore there would be no glitch activity. However, if in- 
dcecd SGR0526-66 is a CFL star, the late formation and a 
resulting recent boost of the transverse velocity would not 
explain its current offset from the nearby young SNR N49: 
Even if the neutrino emission had been anisotropic, the 
SGR can not have moved very far during a few decades, 
and so the late formation scenario would suggest that 
there is no physical association between SNR N49 and 
SGR0526-66. In any case, there is currently no compelling 
observational evidence for such an association (Gaensler 
2002). 

It is important to notice that the CFL phase transition 
will occur only once, so the model discussed in this pa- 
per predicts that there will be only one giant burst. Also, 
note that the probability of observing such an event is 
small, since the mean time between supernovae is about 
100 years and formation of a quark star is maybe even 
rarer. This should be compared to the short time we have 
had gamma-ray observatories. 

It has been demonstrated that CFL stars with very 
short rotation periods should not exist, because they are 
prone to so-called r-mode instabilities (Madscn 2000). 
Since the CFL star we consider has a rotation period of 
the order of a few seconds, this should not be a problem 
for our model. 

Finally, the era of decaying vortices can be estimated 
to last 10 3 to 10 5 years in our model. Even if the formation 
of the quark star happens immediately after (or during) 
the supernova explosion, the SGR should remain active 
long enough to account for the observed ~ 10 s periods. 
A fast spin down occurs if the magnetic field is initially 
large as our model assumes. 

9. Conclusion 

We presented an alternative model for explaining SGRs 
where a newly born quark star experiences bursting ac- 
tivity as it cools below the critical temperature for the 
onset of CFL. In our model, the subsequent magnetic field 
expulsion and reorganization are crucial ingredients in ex- 
plaining the giant bursts and the subsequent weaker ones; 
magnetic field strength and period alone cannot be re- 
sponsible for the unusual properties of SGRs according 
to our model. We also discuss the relation between AXPs 
and SGRs which in our model are separate objects. 
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